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We study the disorder effect of resonant spin Hall effect in a two-dimension electron system with 
Rashba coupling in the presence of a tilted magnetic field. The competition between the Rashba 
coupling and the Zeeman coupling leads to the energy crossing of the Landau levels, which gives 
rise to the resonant spin Hall effect. Utilizing the Streda's formula within the self-consistent Born 
approximation, we find that the impurity scattering broadens the energy levels, and the resonant 
spin Hall conductance exhibits a double peak around the resonant point, which is recovered in an 
applied titled magnetic field. 
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I. INTRODUCTION 

Spin-orbit couplings open a route to control quantum 
electron spin by electric means. One of the efficient meth- 
ods to inject or generate electron spin in non-magnetic 
semiconductors is the spin Hall effect, in which an elec- 
tric current or an electric field may induce a transverse 
spin current in the systems with strong spin-orbit cou- 
plings. Early theories^— proposed that the spin current is 
caused by asymmetric scattering of electrons with spin- 
up and -down in impurity potentials, named as extrinsic 
spin Hall effect. In recent years it was demonstrated that 
the spin-orbit coupling in the electron bands can also 
lead to an intrinsic spin Hall effect in either p-doped or 
n-doped semiconductors.^ Both extrinsic and intrinsic 
spin Hall effects were confirmed experimentally in various 
systems.— 

A two-dimensional electron gas (2DEG) with a Rashba 
coupling was proposed to exhibit an intrinsic spin 
Hall effect 4^ The spin-orbit coupling in 2DEG mod- 
ifies the electron band structure, and may lead to 
interesting magnetotransport properties, such as the 
beating phenomenon in the Shubnikov-de Haas (SdH) 
oscillation. 10 ' 11 When the system is subjected to an exter- 
nal magnetic field, the Zeeman splitting will also change 
the spin-dependent electron bands. The interplay of the 
spin-orbit coupling and the Zeeman coupling produces 
the crossing of electron energy levels. Based on this prop- 
erty, it was proposed that a tiny electric field may remove 
the additional degeneracy of energy levels and produces 
a finite spin current if the Fermi surface sweeps across 
the crossing point of energy levels. As a result, there ex- 
hibits a divergent spin Hall conductance^i^is This 
resonant spin Hall effect was also discussed in p-doped 
systems in a magnetic field . 16 ' 17 

However, impurities in the system make the issue more 
subtle . 18 ' 19 The vertex correction in the self energy turns 
out to cancel the spin Hall conductance even in a weak 
disorder limit in 2DEG with linear Rashba coupling, 
while the spin Hall conductance survives in p-doped 
Luttinger model and the systems with cubic spin-orbit 
couplings. 20 ' 21 ' 22 ' 23 ' 24 ' 25 The disorder effect strongly de- 



pends on the symmetry of the spin-orbital coupling and 
the dispersion. Now whether the resonant spin Hall ef- 
fect can survive in a finite density of impurities becomes 
an issue to be answered. This is the motivation of the 
present work. 

Here we present a full investigation on the disorder ef- 
fect of resonant spin Hall effect in 2DEG with the Rashba 
coupling in a tilted magnetic field. The impurity ef- 
fect is considered by the self-consistent Born approxi- 
mation (SCBA) and the vertex correction in the lad- 
der approximation. We found the impurity effect will 
suppress the resonant spin Hall conductance at the res- 
onant point, and produce a double peak structure of the 
spin Hall conductance around the point. A tilted mag- 
netic field is applied to enhance the effective Zeeman 
splitting ) 26 ' 27 ' 28 ' 29 ' 30 and to recover the effect when the 
energy level splitting is larger than the energy broadening 
by impurities 



II. GENERAL FORMALISM 

A. 2DEG in a tilted field 

We consider a 2DEG in the x-y plane with the Rashba 
spin-orbit interaction in a tilted magnetic field. The per- 
pendicular component of the tilted field is —B±, and the 
in-plane component is chosen to be along the x-direction 
B± tan 9, where 9 is the angle between the field and the 
z-direction. We take the Laudau gauge for the vector 
potential of the field ~B = (B± tan 0,0, -B±). The total 
Hamiltonian including the Zeeman energy is given by 
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where p = —ihsj, m,—e,g s are the electron's effective 
mass, charge, and Lande g factor, respectively, /is is 
the Bohr magneton, A is the strength of Rashba spin- 
orbit coupling, and a are the Pauli matrices. We take a 
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periodic boundary condition along the x direction, hence 
the momentum p x = hk is a good quantum number. 

An analytical solution can be obtained in the case of 
9 = O i 12 i 31 i 32 There were some studies on spin trans- 
port based on the solution i 13 i 33 i 34 Inclusion of the tilted 
field makes the problem much more complicated, and an 
analytical solution is not available at present. In the fol- 
lowing approach, we choose the energy eigenstates for the 
system without the spin-orbit coupling (A = 0) and 9 = 
as a set of basis, 
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where the magnetic length = y / H/eB±, the spin index 
a =1,1 and L x r y \ is the length of the 2DEG, 0„(y) is 
the eigenstate of the n th energy level of a linear oscil- 
lator with the frequency co = eB±/m^ and \a) is the 
eigenstate of spin <j z . 

When the tilt angle 9 = 0, the system can be solved 
exactly . 12 i 31 The eigenvalues of Hq are given by 
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FIG. 1: (Color online) Energy levels as functions of the mag- 
netic field when the tilt angle = 0, the arrow denotes a 
level-crossing point, which develops into a gap when the tilt 
angle increases. The inset shows the gap as a function of the 
tilt angle, reflects the accordance of the numerical calculation 
and the analytic expression Eq. ([7]). The energy has been 
scaled by hui = heB/m. 



where r\ = Xml^jt? and g = g s m/2m e with m e the mass 
of a free electron, s = 1 for n = 0, and s = ±1 for n ^ 1. 
The states & n ks have a degeneracy = L x L y eB jh, 
corresponding to N$ values of k. The eigenstate has the 
form, 
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where f?oi 



= 0, 



and for n ^ 1, 9 ns = arctan(— u n + 
u 2 ) with u n = (1 — g)/\Z8nrj. One of the fea- 



tures of the solution is the crossing of the energy levels 
as functions of the magnetic field, which is caused by the 
competition between the spin-orbit coupling and Zeeman 
energy splitting. For the two levels e„i and e„+i._i, the 
condition for the crossing is determined byi2 



V(l - gf + 8m/ 2 + y/(l - g) 2 + 8(n + l)r? 2 = 2. (5) 

This point is called the resonant point for resonant spin 
Hall effect. 

This additional degeneracy due to the competition be- 
tween the spin-orbit coupling and the Zeeman energy of 
the perpendicular field can be removed by a tilted field. 
In the case of 9 7^ 0, the energy levels can be calculated 
numerically. Using the expression in Eq. (4), we may 
make a truncation approximation by keeping the Lan- 
dau levels with n < N such that the dimensionality of 
matrix is reduced to 2N x 2N. Numerical diagonalization 
of the matrix can give us the energy eigenvalues. 

Alternatively, the gap can also be calculated approxi- 
mately by the degenerate perturbation theory. We take 
the partial Hamiltonian H' = g s [iBB±t&ii9a x /2 as a 
perturbation, and express it in the subspace spanned by 
the two states 1 and $ n +i ^ _i near the resonant 



iA/2 
-iA/2 



where the gap A is 



A = g s ^BB_\_ tan 9 cos 9 n \ sin 9 



n+l, 
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In Fig. [1] we present the energy levels of 9 = as a func- 
tion of B±. The parameters used are A = 9 x l0~ 12 eVm, 
g s = 4, and m = 0.05m e J^- The red arrow denotes a level 
crossing at Bo ~ 2.4T. The insert shows the energy gap 
as a function of the tilt angle 9 with £>_L = Bq. We no- 
tice that the numerical and analytical results are in good 
agreements. 



B. Self-Consistent Born Approximation 

In this section, we briefly review the general formal- 
ism of linear response theory of self-consistent Born ap- 
proximation (SCBA) for electron transport. We shall use 
this technique to investigate the transport properties of 
2DEG with a Rashba coupling in a tilted magnetic field, 
especially near the resonant point. The effect of impuri- 
ties will be taken into account in this formalism. 

We consider a random configuration of impurities with 
short-range potentials V(r) = VS(r — Rj), where 

Rj is the position of the j th impurity. The density of the 
impurities is n, = Ni/(L x L y ). Generally speaking, the 
Green's functions for a specific configuration of random 
potential can be written as Q ± {E) = [E - H - V(r) ± 
z0 + ] _1 , where + and — correspond to the retarded and 
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advanced Green's function, respectively. All transport 
quantities can be expressed in terms of the Green's func- 
tion after averaging all possible configurations of the im- 
purities. Using the conventional perturbation expansion 
with respect to V"(r), we can obtain the Dyson equation 
for the averaged Green's function G^. The impurity ef- 
fect is absorbed by a self-energy function E 1 * 1 as follows, 

G ± (E)^(g ± (E)) c = [E-H -^(E)}, (8) 

where (• ■ ■ } c means the average over all the impurity con- 
figurations. In the SCBA, the self-energy operator can be 
expressed by E± = (VG± (E)V) C M£L2&22. In the repre- 
sentation of the Landau levels, G, £, and V are expressed 
as matrices. For such a spin-independent impurity po- 
tential, previous work o 36 ' 40 proved that the self-energy is 
independent of n and k for a spin-independent Landau 
system. We find that the self-energies for a spin-orbit 
coupling system are independent of n and k 7 
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Here we dropped the index k in G because the averaged 
Green's functions are k-independent. 



C. Streda's Formula for spin Hall conductivity 

With the averaged Green's function in mind, we can 
use the Kubo formula to calculate the linear response of 
any physical quantity O to an external electric field E ext 
in the v direction, 



= lim (O) /E, 



(10) 



As the single-particle version of the Kubo formula, the 
Streda's formula is a conventional and powerful tool to 
study the transport property of 2DEG system under a 
magnetic fieldi 4 ^ At the zero temperature, the formula is 
given by 

o°{E s ) = ieh I dE ■ Tr[O dS *^ v v A{E) 



dE 



-OA(E)v, 



dg-(E). 
dE ■ 



(11) 



where A = (Q~ — Q + )/(2m) is the spectral function and 
v v = -rs[r v ,H] is the velocity operator,. Because there 
are products of two Green's functions in the impurity 
average (• • • ) , the vertex correction has to be included. 
For a specific density of charge carriers, the Fermi energy 
as a function of the magnetic field is determined by 



rEf 

/ dE ■ Tr[A] 
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Usually the Streda's formula is applied to calculate the 
electric conductance by replacing O by an electric current 
operator, J„ = — ev v . In the present work, we intend to 
explore the spin transport in the system. The spin cur- 
rent is defined as j" = (h/A){Vfj,, cr a }, which is a tensor 
determined by both the motion direction of an electron 
and its polarization. In the framework of linear response 
theory, the spin Hall conductivity cr^ v , the ratio of the 
spin Hall current to an external field, can be calculated 
by substituting O = j" (fj, ^ v) in Eq. (fTTj) . The spin 
Hall conductivity comes from the contribution of all the 
electrons below the Fermi level. Opposite to those for the 
conductivity and Hall conductivity, it cannot be reduced 
to a Fermi edge quantity 3 ^ because the spin current is not 
a commutator of any operator and the Hamiltonian.— 
For the purpose of our numerical calculation, we trans- 
form the Streda's formula into the following form, 

= fjj P dE ■ Tr[3»{Kt- K-~% 

(12) 

where 
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K° a are determined in a set of Bethe-Salpeter-like equa- 
tions, 



A-,t- = ^[G+]- 1 F+-+F+- [G-]- 



K++ = 



dE 

+ G+(VK+-V) c G-, 
dG 



dE 



(15) 



dE 



[G+] L F++ + G + {VK++V) c G + . (16) 



where F° a = (g a v v g a j (a, a' G {+, -}) is the vertex 

operator, which satisfy the Bethe-Salpeter equation in 
the ladder approximation, 



T ' =G a (v u + (vF™'v) )G a ' . 



(17) 



These equations can be solved self-consistently and Eq. 
(fl~5|) and (fl"7)) have multi-solutions. From the continuity 
equation for charge current in equilibrium, we can derive 
an auxiliary equation, 

Tr[F+- (£- - £+)] = 0. (18) 

Differentiating Eq. (fT5|) with respect to E leads to an- 
other auxiliary equation, 

Tr[K+~ (E- - £+)] = -Tr[F+ 



,d(E" 



dE 



(19) 



dG a /dE is determined by an another self-consistent 
equation, 

dC a I dC a \ 

— — = —G a G' y + G a ( V —j- — V / G a . (20) 
dE \ dE / „ 
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FIG. 2: (Color online) The spin Hall conductivity as a func- 
tion of the magnetic field around the level-crossing point when 
9 — 0, for various scattering strength F. 



Since the density of states for each Landau level in SCBA 
has a semi-elliptic foriri ) 43 i 44 ' 45 it approaches zero and 
dG a /dE becomes infinity at the edge of each level. We 
find the integrand in Eq. (JT2J) is always convergent in nu- 
merical calculation, because the concurrence of dQ a /dE 
and g a in Eq. JT3]), flU). 



III. NUMERICAL RESULTS 

Now we are ready to calculate the spin Hall conduc- 
tivities numerically. In this paper, the electron density is 
fixed at n e = 2.9 x 10 15 ?ti~ 2 . This value of n e promises 
that the Fermi level is located near the resonant point 
with the filling factor v = n e /{N^/ L x L y ) = 5 while the 
magnetic field sweeps over the point, i.e., B± = Bq as in- 
dicated in Fig. [TJ The other parameters are as the same 
as those used in Fig. [TJ The density of impurities n.j and 
the strength of the impurity potential V are combined in 
one parameter "scattering strength" T = niV 2 m/(2Trh 2 ). 
We assign T various values to investigate the impurity ef- 
fect. 



A. Disorder effect of Resonant spin Hall effect 

We first discuss the disorder effect of spin Hall conduc- 
tance, especially near the crossing point. We apply the 
formula in Eq. lfT!?)) to calculate the spin Hall conductiv- 
ity a z xy around the resonant point Bq = 2AT numerically 
for various strengths of disorder. Numerical results are 
plotted in Fig. [51 The dashed curve for V = is from 
the solution in Ref. [l2j|. The key feature of the disor- 
der effect is the suppression of the spin Hall conductivity 
at the resonant point. The large spin Hall conductance 
exhibits when the field deviates from the resonant point 
and forms a double peak structure. The weight of the 
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FIG. 3: (Color online) The distribution of the spin Hall con- 
ductivity on the electron's energy, da% y /dE, as a function of 
the energy E and the magnetic field B± when the tilt angle 
9 — 0, with the impurity strength F = 1/32/j.eV . The energy 
E has been scaled by huj. 



spin Hall conductivity increases as the impurity strength 
decreases, which reflects the intrinsic properties of the 
resonance. 

To understand the suppression of resonant spin Hall 
conductance, we plotted in Fig. [3] the distribution of the 
spin Hall conductivity on the electron's energy, da xy /dE, 
around the crossing point. The total spin Hall conductiv- 
ity <j xy {Ef) = J^£ c (da* y /dE)dE are contributed by all 
the electron states under the Fermi level. It was observed 
that the energy levels are broadened due to the impurity 
scattering, and the distribution of the spin Hall conduc- 
tivity is inhomogeneous. So the magnitude and even the 
sign of the spin Hall conductivity can be varied by the 
Fermi level or the electron density. In the clean limit, a 
tiny external electric field can open a gap between the 
crossing levels, which leads to spin Hall conductance di- 
vergent. However, after the impurity scattering is taken 
into account, a tiny external field cannot open an energy 
gap any more because of the level broadening. As a re- 
sult, the spin Hall conductance in a weak field limit will 
be suppressed. However, once the external field becomes 
stronger than the level broadening, a large spin Hall con- 
ductance will appear. This can be seen from the case 
that the magnetic field deviates from the crossing point, 
i.e., the additional degeneracy of the two levels will be 
lifted, and a strong spin Hall conductance recovers. This 
is the physical origin of the double peak structure of the 
resonant spin Hall conductance. It is worth stessing that 
this suppression of resonant spin Hall conductance is dif- 
ferent from the case in the absence of the Zeeman term. 
The Zeeman splitting may produce a non-zero spin Hall 
conductance in the Rashba system^ 
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B. Effect of a tilted field 

To further illustrate the formation of the resonant spin 
Hall effect, we investigate the effect of the tilted magnetic 
field. Fig. 2] shows the dependence of the spin Hall con- 
ductivity on the tilt angle near the resonant point. For 
the purpose of numerical calculation, we take the scatter- 
ing strength T = 1/lQ/ieV. As the tilted angle increases, 
the spin Hall conductivity at the resonant point will in- 
crease very quickly, and the two peaks finally integrate 
into one. After that point, the spin Hall conductivity 
begins to decrease. These behaviors can be understood 
as the competition between the disorder broadening of 
the energy levels and the degeneracy lifting by the tilted 
field. 

In the clean limit, the tilted field will remove the de- 
generacy of the energy crossing levels as shown in Fig. 
[TJ We can estimate the peak height of the spin Hall 
conductivity as a function of the the tilt angle by a per- 
turbation calculation adopted in Sec. Ill Al Diagonalizing 
the truncated two- level Hamiltonian in Eq. ©, we get 
the modified eigenstates ty± = ($ n fci ± i&n+i.k.-i) /v2, 
and the energy correction E± = ±A/2 with Eq. (JT)). If 
the Fermi level just lies between the energy levels, the 
spin Hall conductivity is mainly attributed to which 
can be calculated by the Kubo formula ) 12 ' 13 



<4 1} obo) = E( 
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\{n + 1) cos 2 6 nl cos 2 9 n 
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he 2 B 
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(21) 



where y = (mv x \\ = o — Hk)/eB. 

The impurity scattering will cause the level broaden- 
ing. If the level broadening is larger than the gap caused 
by the tilted field, the impurity effect is dominant. Oth- 
erwise the tilted field effect will be dominant. The level 
broadening is characterized by the half- width of the semi- 
elliptic density of states of the Landau levels, which can 
be estimated approximately by^ 7 = 2^mV 2 N$>/ L x L y . 
When A > 27, the effect of the tilted field turns out to be 
dominant. The insert of Fig. [4] presents the data of the 
peak height of the spin Hall conductivity as a function 
of 9 from the numerical calculation (diamond) and from 
the analytic formula (solid line) in Eq. ([2"!]) . They are in 
a good agreement when 9 10°, which is very close to 
the estimated value 9 ss 8° from A = 27. 



2DEG system with the Rashba interaction in a tilted 
magnetic field. Considering the vertex corrections in the 
self energy, we find that the main effect of the impu- 
rity scattering is to broaden the Landau levels. In the 
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FIG. 4: (Color online) Spin Hall conductivity as a function of 
the perpendicular component of the magnetic field for various 
tilt angles 9. The insert compares the peak height of the 
spin Hall conductivity as function of 6 from the numerical 
calculation (diamond) and from the analytical formula Eq. 
(pTj) (solid line). The scattering strength V = l/16fj,eV. 



framework of linear response, the electric field is taken 
to approach zero, and the energy splitting caused by the 
electric field is always less than the broadening of the 
Landau levels. Thus a tiny external field cannot remove 
the additional degeneracy of the energy levels at the res- 
onant point. As a result, the spin Hall conductance will 
be suppressed at the point. When the magnetic field 
slightly deviates the resonant point or a tilted field is ap- 
plied, the degeneracy will be removed, a large spin Hall 
conductance will be recovered. The spin Hall conduc- 
tance exhibits a double peak around the resonant point. 
From the effect of a tilted field, we believe that a finite 
electric field, if it is strong enough to overcome the energy 
level broadening, will recover the spin Hall effect even at 
the resonant point. This is quite different from the dis- 
order effect of spin Hall effect in the Rashba system in 
the absence of magnetic field.— 
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